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ABSTRACT 

A  mathematical  model  of  a  vehicle  being  brought  to  rest  in 
such  a  manner  as  to  cause  a  skid  is  devised  and  simulated  on  an 
analog  computer.   This  model  includes  the  effects  of  the  brake 
line  pressure  upon  the  braking  of  the  vehicle  and  also  the  effects 
of  the  coefficients  of  friction  between  the  tire  and  road  surface 
and  between  the  brake  shoe  and  brake  drum.   Several  controls  to 
correct  for  the  occurrence  of  a  skid  are  simulated  on  the  analog 
computer  and  a  detailed  analysis  of  the  braking  system  under  the 
influence  of  each  control  is  performed.   Lastly,  a  scheme  for 
implementing  the  control  which  minimizes  the  stop  time  and  the 
degree  of  skid  is  presented. 
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1.   INTRODUCTION 

Skidding  of  a  vehicle  due  to  an  excessive  application  of 
brake  force  for  the  road  conditions  prevailing  causes  many 
thousands  of  accidents  every  yeara   Rapid  and  heavy  application 
of  brakes  in  a  panic  stop  situation  often  causes  wheel  lock  up 
and  subsequent  loss  of  control  of  the  vehicle.   This  also  occurs 
in  wet  or  icy  road  conditions. 

The  minimum  stopping  distance  possible  for  a  vehicle  is 
determined  by  the  coefficient  of  friction  between  the  tire  and 
the  road  surface.   When  a  vehicle  is  just  on  the  verge  of  a  skid 
the  frictional  force  between  the  tire  and  road  surface  is  at  its 
greatest  value,  and  consequently  the  deceleration  of  the  vehicle 
is  at  a  maximum,  i.e.,  equal  to  |i  g,  where  p.   is  the  coefficient 
of  static  friction  between  the  tire  and  the  road  surface  and  g 
is  the  deceleration  of  gravity  (see  Appendix  A  for  a  proof  of 
this  statement).  When  skidding  occurs  the  tire  slides  along  the 
road  surface.   This  causes  a  reduction  in  the  coefficient  of 
friction  and  a  corresponding  reduction  in  the  frictional  force 
between  the  tire  and  road  surface,  thus  decreasing  the  magnitude 
of  the  vehicular  deceleration.   The  optimum  braking  situation 
would  thus  be  attained  if  the  vehicle  could  be  kept  just  at  the 
verge  of  a  skid.   If  a  skid  could  be  detected  exactly  at  the 
instant  of  occurrence  and  a  compensating  force  were  automatically 
applied  to  eliminate  the  skid  before  wheel  lock  up  occurred,  then 
nearly  ideal  braking  conditions  could  be  realized  under  all  road 
conditions . 


Attempts  have  been  made  previously  to  develop  "anti-skid" 
devices,  and  some  such  devices  are  presently  in  use  on  certain 
aircraft  and  many  trucks.   Previous  efforts  in  this  direction 
include  the  Dunlop  System  of  England,  the  Maxaret  Unit,  the  Anti 
Skid  Device,  the  Westinghouse  Decelostat  Controller,  the  Hydro- 
Aire  Hytrol  system,  and  the  Anti-Skid  Braking  system  of  Russia's 
Tu  -  104  Jet  Liner. [l]   These  devices  all  work  satisfactorily  to 
a  greater  or  lesser  extent,  but  most  have  the  drawback  that 
operation  depends  upon  the  sensing  of  either  a  complete  wheel 
lock  up  or  the  attainment  of  a  pre-determined  rate  of  wheel 
deceleration,  at  which  time  a  bang-bang  type  of  control  is  initi- 
ated. Additionally,  many  of  these  devices  require  expensive, 
heavy  installations,  rendering  them  largely  unsuitable  for  eco- 
nomical installation  and  operation  on  commercial  vehicles. 

This  thesis  presents  an  analysis  and  evaluation  of  several 
electrically  actuated  alternate  proposals  designed  to  automatically 
compensate  for  an  overbraking  situation  in  order  to  minimize  the 
degree  of  skidding  and  more  nearly  optimize  the  braking  of  a  vehicle, 
A  mathematical  model  for  the  system  was  developed  and  simulated  on 
an  analog  computer,  and  three  types  of  skid  correcting  controls 
were  devised  and  tested.   One  form  of  control  which  linearly 
decreases  the  net  force  applied  to  the  brake  shoe  until  the  skid 
is  eliminated  and  then  permits  a  linear  increase  of  the  net  braking 
force  until  a  skid  once  again  occurs  was  tested.   This  is  called 
the  double  ramp  control.   Two  forms  of  bang-bang  controls  were  also 
tested.   Finally  a  triple  ramp  control  which  continually  "samples" 
the  road  surface  to  determine  the  optimum  brake  force  was  devised 
and  tested. 
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Complete  data  -were  obtained  for  each  type  of  control  and  a 
detailed  analysis  was  made  of  the  behavior  of  the  braking  system 
under  the  influence  of  each  control. 

Lastly,  a  scheme  for  implementing  the  triple  ramp  control 
(which  proved  to  be  superior  to  any  others  devised)  is  advanced, 
although  the  physical  system  was  not  actually  built. 
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2.   DESCRIPTION  OF  THE  PROBLEM 

For  the  purposes  of  this  paper,  a  skid  will  be  defined  as 
that  condition  which  exists  when  the  actual  linear  velocity  of 
the  frame  of  the  vehicle  in  question  is  greater  than  the  product 
of  the  angular  velocity  of  the  wheel  and  the  radius  of  the  wheel, 
In  terms  of  the  symbology  used  throughout  this  paper 


e  =  V  -  R-cO,  (1) 


e>0     when  v>R.cJ  (2) 


and  e   =  0     when  v  =  R,^,  (3) 


where  v  represents  the  linear  velocity  of  the  frame,  R..  represents 
the  radius  of  the  wheel  and  CO  is  the  angular  velocity  of  the  wheel. 
The  term  e  corresponds  to  an  error  signal,  the  magnitude  of  which 
is  a  measure  of  the  degree  of  skidding  present.   The  product  R  cO 
will  be  referred  to  as  the  velocity  measured  by  the  rotation  of  the 
wheel.  When  e  =  0  there  is  no  skidding  and  when  e  ^>  0  the  vehicle 
is  moving  faster  than  a  point  on  the  rim  of  the  wheel,  i.e.,  the 
wheel  is  tending  to  lock  up  and  a  skid  is  occurring.   The  condition 
when  R..cJ>V  would  occur  when  the  wheels  are  spinning,  as,  for  example, 
when  trying  to  accelerate  on  an  icy  road. 

A  proposed  means  of  sensing  the  linear  velocity  of  the  vehicular 
frame  is  to  integrate  the  output  from  a  linear  accelerometer  attached 
to  the  axle  housing  of  the  vehicle  (to  minimize  effects  of  vehicle 
tilt  during  braking),  and  to  compare  the  integrator  output  with  the 
velocity  measured  by  the  rotation  of  the  wheel.   The  initial  condition 
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input  to  the  integrator  would  be  the  velocity  measured  by  the 

rotation  of  the  wheel  at  the  instant  of  application  of  the  brake 

force. 

When  an  error  signal  occurs,  indicating  skidding  is  present, 

an  automatic  corrective  force  (F  )  would  be  generated  to  reduce 

the  total  force  F  applied  to  the  brakes.   F.   represents  the  force 

rr  in   r 

applied  by  the  vehicle  operator  to  the  brakes  and  is  considered 
constant  for  all  subsequent  analyses. 

These  concepts  are  illustrated  in  Fig.  1  where  a  block  diagram 


of  the  proposed  system  is  shown. 
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Fig.  1.   Block  diagram  of  the  proposed  braking  system. 


From  elementary  physics  it  is  known  that  the  coefficient  of 
friction  between  two  surfaces  is  dependent  upon  the  relative  velocity 
between  the  surfaces.   A  sketch  of  the  behavior  of  the  coefficient 
of  friction  as  a  function  of  the  relative  velocity  between  the  two 


surfaces  is  shown  in  Fig.  2. 
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Fig.  2.   Coefficient  of  friction  as  a  function  of  relative  velocity. 

Just  before  relative  motion  begins,  the  coefficient  of  friction 
between  two  surfaces  is  at  its  maximum  value.   Prior  to  reaching 
this  maximum  value,  the  coefficient  of  friction  is  a  variable  that 
adjusts  itself  in  such  a  manner  so  that  the  force  tending  to  cause 
a  relative  motion  is  just  compensated  by  the  reactive  force  of 
friction  (this  is  before  any  relative  motion  takes  place) .   After 
sliding  motion  starts  (as  between  a  brake  shoe  and  drum,  or  between 
a  tire  and  road  surface  for  a  skidding  situation)  the  coefficient  of 
friction  |i  varies  as  shown. 

A  piecewise  linear  approximation  to  this  behavior  is  shown  in 
Fig.  3. 


k 


/*• 


if^/^s 


/u-s>A'o 


vc 


RELATIVE- 
VELOCITY 


Fig.  3. 


Piecewise  linear  approximation  to  the  behavior 
of  the  coefficient  of  friction. 
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The  term  u.  corresponds  to  the  coefficient  of  static  friction 

s 

which  is  a  maximum  just  before  relative  motion  begins,  and  which  is 
always  greater  than  the  coefficient  of  kinetic  or  sliding  friction. 
The  equation  for  this  piecewise  linear  approximation  is 


jU  =  /J-o    ,  -for        0<V<£Vc 


(4) 


/<=/*• -K,'(V-VC)  }+or     Vc^  V^(vc   -hSf- j    . 


(5) 


The  above  relationship  will  be  used  to  represent  the  coefficient 

of  friction  between  the  brake  drum  and  brake  shoe  of  the  moving  vehicle. 

The  break  point,  v  ,  for  the  coefficient  of  friction  between  the  tire 

c 

and  the  road  surface  will  be  assumed  to  be  high  enough  so  that  for 

normal  speeds  the  coefficient  of  kinetic  friction  is  some  constant 

|i  .   As  will  be  seen,  this  is  not  a  restrictive  assumption,  since  the 
K. 

purpose  of  the  corrective  forces  to  be  applied  is  to  minimize  the  degree 
of  skid  experienced.   Thus  the  relative  velocity  between  the  tire  and 
road  surface  would  be  small  enough  to  warrant  considering  the  coefficient 
of  kinetic  friction  to  be  a  constant. 

Consider  the  idealized  brake  drum  and  shoe  shown  in  Fig.  4. 
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Fig.  4.   Diagram  of  brake  drum  and  shoe. 
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In  Fig.  4  the  wheel  is  assumed  to  be  rotating  in  a  clock- 
wise direction  with  angular  velocity  uu  and  angular  acceleration 
oo  with  the  positive  sense  as  shown.   F  represents  the  total  (net) 

force  acting  on  the  brake  shoe,  and  is  the  difference  between  F. 

m 

(applied  by  the  brake  pedal)  and  F   (the  corrective  force).   F 
is  the  force  of  friction,  the  value  of  which  depends  upon  whether 
a  skid  is  present.   The  equations  describing  this  system  are  pre- 
sented below. 

F  =  F.      -   F  (6) 

in  c 

Net   torque   =  T   =  R0   Fp.   -  R.F_  (7) 

Ff  =      f~mv     ;    for    e  =  O  (8) 

I    /*k  fn  g   ,  for     e>  O  (9) 

where  m  =  the  weight  of  the  vehicle  acting  through  the  axis  of 

2 

the  wheel  divided  by  the  acceleration  of  gravity  g  (=  16  ft. /sec.  ), 

,,  _    ( M-°-   K/R*  (^-^c)j  for     CO^Cdc  (10) 

I  /j.0      ;    -for      Co  ^  C0o  (u) 

i.e.,    |X   is   the   coefficient   of    friction  between   the  brake   drum  and 
brake   shoe    (Fig.    3).      Also,    \i,      is   the   coefficient   of   friction  between 
the   tire  and  the  road  when  a   skid   is   present. 
Let     K,    R^  =    K,  ,    so   that   Eq.    (10)    becomes 

=  f /*<>-  k,  (oo-co^)  )  f  or     o;>60c  (12) 

[  /U.0    ,     -Pot       CO  <C  U)c  (13) 

and  Eq.    (7)   becomes 

T=  R£  F/u-   -R,  Ff  =  -  Jco  (14) 

where  J  =  the  moment  of  inertia  of  the  wheel. 
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Consider  the  case  where  e  =  0,  i.e.,  skidding  is  absent,  so 

F+  =    -   hi  V  ^  /Us  J  J  (15) 

with  the  equality  holding  when  the  system  is  just  at  the  verge  of 
a  skid.   The  equation  of  motion,  (14),  becomes 

-Jo>  =  R,  fn  V  +  R£  F/^  .  (16) 

Since  e  =  0,  then 

v  =   R{  cb    .  (17) 

Substituting    (17)    into    (16)   yields 

-J co  =  R,£  m  cb  +  Ra  F/M  .  (18) 

Solving    (18)    for  <J0  yields 
J+  RfVn 


(20) 


When  e  >  0,    i.e.,    skidding  is   present,    and 

Ff  =  ^  k  ^  J     j    whei-e     //~*</u-s    > 
then   (14)   becomes 

-J,  lc  =  Re  F'/c  -  R  ,  /6K  rn  g    ;  (21) 


so 

6j  =  R/*n  ^3_  -    R&F/y  .  (22) 

Note   that    for   each   of  the  relations   for  <J0   the  value   of  y.   is 


a    function   of  U) .      In  the   case  where  uo  <  ou    ,    then  p.   =  p,    ,   and   the 

c  o 

equations   above  apply  with  \i  replaced  by  p.    .     When  UJ  ^  w    ,    then 

/^/<o    -    K,  (a)-u)c)  .  (23) 

The   complete   set   of  equations    is   indicated  below: 
When  e  =  0 ,    co  <  CO  c   ; 

&  =      ~  R*    F        yU0     .  (24) 
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When  e  >  0,    CO  <  CO  c  I 

cb  m    R./^Krng     -    RaEaao    .  (25) 

7  7 

When  e  =   0,   CO  ^  Cdc  I 


Cb  =    R^FK,     ^    _    Rg  FK£     ,  (26) 

7>Rfm  jtRjam 

When  e  >   0,  CO>  OdQ  '. 

co  =    Ra  FK,    ^    -     RaFKe     +■     R    AKh^q     .  (27) 

J  J  J 

In  the  last  two  relations, 

The  analog  computer  simulation  of  this  system  of  equations  is 
shown  in  Fig.  5.   Note  that  Appendix  B  contains  a  brief  description 
of  the  analog  computer  symbols  used  in  this  paper  and  an  explanation 

of  the  logic  elements  used  to  control  various  portions  of  the  simu- 

,   .    C2,3,¥,5J 
lation. 

To  determine  the  equations  governing  the  forward  motion  of  the 

vehicle,  recall  that  the  maximum  possible  deceleration  for  a  given 

coefficient  of  friction  p, '  between  the  tire  and  road  surface  is 

-xV«x  =Xg    J  (28) 

where  Ll     -    JUL  s  or     M.  ^  . 

Thus,  when  e  =  0,  i.e.,  no  skidding, 

-Xs-R.i^/tsJ  j  (29) 

where  M-  =  the  coefficient  of  static  friction  between  the  tire  and 

s 

the  road  surface.   Also,  before  a  skid  occurs, 

X  =  R.  CO    =  V  .  (30) 
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After  a  skid  occurs, 


(31) 


where  M-   =  the  coefficient  of  kinetic  friction  between  the  tire 
K 

and  road  surface  which  is  considered  constant  for  this  simulation, 

and   /-<  K  <  /"-  S  • 
Thus,  when  e  >  0, 


V  =  -yU  K  3  t   +  V0    ) 


where  V  is  the  velocity  of  the  vehicle  when  the  skid  begins. 

The  analog  computer  implementation  of  these  equations  and 
ideas  is  presented  in  Fig.  6,  which  provides  a  simulation  for  V 
under  conditions  of  skidding  or  skidless  braking. 
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at     the       ow-t  piLt  . 


Fig.  6.  Analog  computer  simulation  of  vehicular  velocity. 

The  composite  analog  computer  set  up  for  both  the  brake  drum 
and  shoe  and  the  vehicular  velocity  is  shown  in  Fig.  7.   This  is 
intended  to  represent  a  conventional  braking  system  in  existence 
on  present  day  vehicles  and  no  corrective  force  is  included  in 
this  diagram. 
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3.   PROPOSED  SKID  CORRECTING  CONTROLS 

Any  form  of  corrective  control  must  be  initiated  by  the 
occurrence  of  a  skid,  which  in  turn  is  indicated  by  the  presence 
of  a  positive  error  signal.   The  purpose  of  the  corrective  control 
must  be  to  lessen  the  net  force  F  applied  to  the  brake  shoe  when 
a  skid  is  sensed.  When  the  skid  has  ceased,  the  corrective  force 
F  must  decrease  to  permit  a  reapplication  of  the  full  braking 
force  to  the  brake  shoe.   Thus  the  value  of  F  will  depend  upon 
whether  a  skid  is  present  or  not. 

The  first  form  of  control  to  be  tried  is  called  a  double  ramp 
control,  and  a  typical  shape  for  the  control  is  shown  in  Fig.  8. 


e>o 


time 


Fig.  8. 


Double  ramp  control. 


An  examination  of  Fig.  8  indicates  that  the  maximum  value 

attained  by  F  is  F.   (the  force  applied  by  the  vehicle  operator). 

Since  F  =  F.   -  F   ,  and  since  there  can  never  physically  be  a 
in    c  v  j  j 

negative  force  on  the  brake  shoe,  F  ^  F.  .   The  value  of  K_  is 

'cm  3 

not  necessarily  equal  to  K, . 

An  analog  computer  simulation  of  the  double  ramp  control  is 
shown  in  Fig.  9. 
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Fig.  9.  Analog  computer  simulation  of  the  double  ramp  control. 

The  integrator  indicated  in  Fig.  9  may  be  in  the  compute  mode 
for  the  entire  problem  simulation,  due  to  the  presence  of  the  feed- 
back limiter  at  the  output  of  the  integrator  which  prevents  F 
from  ever  going  negative. 

A  bang-bang  control  was  tried  for  the  second  form  of  corrective 
force.   Since  any  physically  realizable  bang-bang  control  would  have 
either  a  dead  zone  or  a  time  delay,  one  control  with  a  dead  zone 
and  one  with  a  time  delay  were  tried.   The  expression  "dead  zone" 
is  meant  to  imply  that  there  is  some  positive  value  of  error  volt- 
age above  which  the  corrective  force  is  full  on  until  the  error  is 
once  again  reduced  to  zero.   The  analog  computer  implementation  of 
the  bang-bang  control  with  a  dead  zone  is  shown  in  Fig.  10,  where  at 
is  the  value  of  the  dead  zone  voltage. 
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Fig.  10.   Analog  computer  implementation  of  the 
bang-bang  control  with  a  dead  zone. 
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The  expression  "time  delay"  as  used  here  implies  that  the 
error  exists  for  some,  time  T1  before  corrective  action  is  initiated. 
In  order  to  simulate  the  time  delay,  a  Wood  fourth-order  approximant  [2] 
time  lag  circuit  was  set  up  to  delay  the  error  voltage  by  one  second 
(real  time  delay) ,  which  in  problem  time  corresponded  to  a  delay  time 
of  0.1  second.   The  analog  computer  diagram  for  the  bang-bang  control 
with  a  real  time  delay,  with  the  potentiometer  settings  needed  to 
produce  a  one  second  real  time  delay,  is  shown  in  Fig.  11. 
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Fig.  11.  Analog  computer  implementation  of  the 

bang-bang  control  with  a  real  time  delay, 


The  third  form  of  control  to  be  simulated  is  called  the  triple 
ramp  control,  and  a  typical  plot  of  corrective  force  F  versus  time 
plus  the  complementary  plot  of  F  (=  F.   -  F  )  versus  time  is  shown 


in  Fig.  12. 
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Fig.  12.   Triple  ramp  control. 

The  triple  ramp  control  is  of  the  following  form.   F  =  0 

until  a  skid  occurs.  When  a  skid  occurs,  e  >  0  and  F  begins 

c   ° 

rising  rapidly  with  a  slope  K~  to  a  maximum  value  of  F.   or  until 

e  =  0.  When  e  =  0,  F  then  begins  decreasing  with  a  slope  equal 

to  -  K,  until  e  >  0  once  more.   Up  to  this  point  the  control  is 

identical  to  the  double  ramp  control.   However,  the  second  time 

that  e  >  0  a  separate  monitoring  device  with  a  "memory"  (a  Mylar 

capacitor,  for  example)  holds  the  voltage  at  the  value  corresponding 

to  the  corrective  force  at  which  the  skid  occurred. 

At  the  instant  when  e  >  0  for  the  second  time  in  Fig.  12,  F 

°    '   c 

would  once  again  increase  rapidly  with  slope  K„  until  e  =  0.   Then 
F  would  again  begin  decreasing  with  slope  -  K  until  it  is  just 
above  the  corrective  force  at  which  skidding  occurred  previously. 
Then  F  would  be  allowed  to  decrease  more  slowly  with  a  slope  -  K 
(thus  permitting  a  gradual  increase  in  F  until  a  skid  once  again 
occurs).  The  cycle  then  repeats  itself. 
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The  design  of  this  particular  type  of  control  was  motivated 

by  the  considerations  expressed  in  the  Introduction,  i.e.,  that 

the  optimum  braking  situation  would  be  realized  if  the  vehicle  could 

be  kept  just  at  the  verge  of  a  skid.   The  storage  of  the  value  of 

F  at  which  a  skid  occurs  endows  the  control  with  a  form  of  "intel- 
c 

ligence"  which  permits  the  net  force  applied  to  the  brake  drum  to 

be  just  less  than  the  optimum  force.   The  gradual  increase  in  F 

permitted  by  the  third  portion  of  the  control  (with  a  slope  of  -  K  ) 

is  necessary  to  permit  the  control  to  adapt  to  improving  (i.e.,  less 

slippery)  road  conditions.   The  first  cycle  of  the  control,  where 

F   increases  rapidly  with  a  slope  of  K„  and  then  decreases  with  a 

slope  of  -  K,  until  the  skid  occurs  once  again,  is  necessary  to 

permit  the  control  to  "test"  the  road  surface  for  the  net  force 

required  to  just  cause  a  skid.   The  analog  computer  implementation 

of  this  control  is  shown  in  Fig.  13. 

Integrator  No.  1  provides  the  ramp  with  a  slope  of  L,  -  K,  , 

or  -  K  as  appropriate.   This  integrator  is  in  the  compute  mode  for 

the  entire  simulation  time  and  F  is  kept  from  going  negative  by  the 

lower  bound  of  the  feedback  limiter.   The  upper  bound  on  the  limiter, 

and  consequently  on  F  ,  is  determined  by  F.  .   Potentiometer  No.  1, 

labeled  G.R./10,  determines  the  value  of  F  at  which  the  slope 

c 

changes  from  -  K,  to  -  K  ;  the  ratio  of  F  /„\/F  ,.v  (see  Fig.  12) 

is  termed  the  Correction  Ratio  (C.R.).   It  is  this  ratio  which 

determines  the  setting  of  potentiometer  No.  1.  When  integrator  No.  2 

is  in  the  track  (reset)  mode  switch  No.  1  is  closed  (FF  =  1)  and  we 

desire  the  slope  of  F  to  be  either  K„  or  -  K.  (but  not  -  K  ),  since 

c  3       4  c 

this  is  the  time  period  between  the  first  and  second  occurrences  of 
the  positive  error  signal.   Thus  to  positively  prevent  relay  No.  1 
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Fig.   13.     Analog  computer  implementation  of  the  triple 
ramp  control. 
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from  throwing  during  this  phase  of  the  cycle  potentiometer  No.  2 
is  set  lower  than  potentiometer  No.  3  to  insure  that  C  =  1  and 
C  =  0,  where  1  and  0  correspond  to  a  logical  true  and  false, 
respectively . 

The  flip  flop  (the  output  of  which  is  denoted  by  FF  and  FF) 
is  an  R  -  S  flip  flop,  and  is  used  to  control  the  mode  of  integrator 
No.  2  as  well  as  the  state  of  switches  No.  1  and  2. 
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4.   THE  UNCONTROLLED  SITUATION 

The  analog  computer  simulation  of  the  uncorrected  braking 

of  the  vehicle  (shown  in  Fig.  7)  was  set  up  on  the  COMCOR  Ci 

5000  Analog  Computer,  using  the  following  arbitrary  parameters: 

K,  =  0.01 

0)  =60  rad./sec. 
c 

[i     =  0.8 
o 

Thus  p.  =  0.8  -  0.01  (U)  -  60),  uo  ^  u> 

and  |o.  =  0.8,  ou  <  u)  , 

c 

where  |i  =  coefficient  of  friction  between  the  brake  drum 

and  shoe. 

la   =  0.75  =>  static  friction  between  tire  and  road  surface. 

s 

10.   =  0.5  ==>  kinetic  friction  between  tire  and  road  surface. 
K 

V  =  75  ft. /sec.   =  initial  vehicular  velocity, 
o 

uo   =75  rad/sec.   ==>  initial  rotational  velocity  of  the  wheel, 
o 

R-  =  1  ft. 

R2  =  .6  ft. 

J  =  25  ft.  lb.  sec.2 


M  =  93.75  slugs 

(Weight  =  mg  =  2000  lbs.  =  weight  of  vehicle  acting 
through  axis  of  wheel). 

A  7000  lb.  input  force  (corresponding  to  F.  )  was  simulated  as 

a  70  volt  step  input  occurring  at  time  =  0  sec. 

Six  different  variables  were  plotted  for  each  run  on  a  Brush 

Instruments  Mark  200  (6  channel)  Recorder.   They  were 

(JJ  ==>  the  angular  acceleration  of  the  wheel,  in  rad/sec.' 

to  =>  the  angular  velocity  of  the  wheel,  in  rad/sec. 


29 


V  ==>  the  actual  linear  velocity  of  the  vehicle, 
in  ft/sec. 

F  =>  the  net  force  applied  to  the  brake  shoe 

(=  F.   -  F  ),  in  lbs. 
in    c 

e  =>  the  error  signal  (=  v  -  R,UJ),  in  ft/sec. 

F  ==>  the  corrective  force,  in  lbs. 
c 

The  uncorrected  system  response  of  the  vehicle,  i.e.,  the 
response  obtained  when  no  corrective  force  is  applied,  is  shown  in 
Fig.  14.   From  examination  of  the  uncorrected  system  plot  it  is 
seen  that  the  wheel  locks  in  0.96  sec.  and  skidding  occurs  0.16 
sec.  after  application  of  the  braking  force.   The  skid  is  present 
until  the  vehicle  comes  to  rest  4.64  sec.  after  the  force  is  initially 
applied. 

Measuring  the  slope  of  the  velocity  output  plot  gives  an 

2 
average  deceleration  of  16  ft/sec.   (excluding  the  first  0.1  sec. 

of  the  problem,  where  skidding  does  not  occur).   This  is  expected, 

2 

since  with  p.   =  0.5  the  deceleration  during  a  skid  is  16  ft/sec. 

K 

(=  0.5  x  32  ft/sec.2). 

To  completely  understand  the  shapes  of  the  various  plots 
displayed  in  Fig.  14,  recall  that  the  coefficient  of  friction  between 
the  brake  drum  and  brake  shoe  is  of  the  general  shape  indicated  in 
Fig.  3,  with  a  break  point  at  uo  =60  rad/sec.   The  initial  rotational 

velocity  of  the  wheel  was  oo   =75  rad/sec,  and  thus  at  the  time  of 

J  o 

simulated  application  of  the  brakes  the  coefficient  of  friction  \i> 
between  the  brake  drum  is 

M.  =  0.8  -  0.01  (75  -  60)  =  0.65. 
The  value  of  \i   will  increase  as  the  relative  velocity  between 
the  brake  shoe  and  brake  drum  decreases,  until  it  has  attained  a 
maximum  value  of  0.8.  After  the  problem  begins  there  is  a  time  lapse 
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Fig.   14.     Uncorrected  braking  situation 
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of  0.16  sec.  before  skidding  starts.  During  this  initial  0.16 
sec.  |i  is  increasing  since  the  vehicle  is  slowing  down,  However, 
p.  is  sufficiently  small  during  this  time  interval  so  that  the  net 
torque  T  (=  R-Fij,  -  R  F  )  acting  on  the  wheel  is  not  sufficient  to 
cause  a  skid.  When  the  velocity  has  decreased  to  70  miles/hr., 
and  id  is  down  to  70  rad/sec.  (since  a  skid  has  not  occurred  yet), 
the  value  for  p.  has  increased  to 

M.  =  0.8  -  0.01  (70  -  60)  =  0.7. 

This  value  for  |i  is  sufficiently  large  to  cause  the  net  torque 

to  be  just  equal  to  that  value  necessary  to  cause  a  skid.   At  the 

instant  a  skid  occurs  the  coefficient  of  friction  between  the  tire 

and  road  surface  drops  abruptly  from  |i  =  0.75  to  |i   =  0.5.   F  also 

S  K  t 

decreases  abruptly  due  to  the  occurrence  of  the  skid  and  the  net 
torque  increases  rapidly,  causing  the  wheel  to  decelerate  rapidly 
and  lock  after  0.96  sec.  total  elapsed  time. 

Note  that  the  plot  for  the  radial  acceleration  of  the  wheel  00 

indicates  that  at  time  =  0  the  acceleration  drops  sharply  to  a 

2 

negative  value  (-  23  rad/sec.  )  and  then  slowly  becomes  more  negative 

until  the  time  when  a  skid  occurs.   During  this  time  the  wheel  is 
still  slowing  down  gradually,  and  the  deceleration  (or  negative 
acceleration)  is  increasing  slowly  as  the  value  of  p.  (between  the 

brake  shoe  and  brake  drum)  increases.  At  the  point  when  the  value 

2 
of  10  becomes  just  less  than  -  24  rad/sec.  ,  the  skid  occurs,  and  to 

2 
abruptly  becomes  much  more  negative  (-  76  rad/sec.  )  and  decreases 

2 
rapidly  until  it  reaches  a  constant  -  94  rad/sec.    The  significant 

factor  to  be  noted  here  is  that  the  skid  occurs  just  when  uo  - 

-  24  rad/sec.  .   This  corresponds  to  a  value  of  V  =  -  24  ft/sec. 

(since  R.  =  1  ft.,  and  a  skid  has  not  occurred  yet).   This  is  in 
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agreement  with  the  previous  statement  that  the  maximum  deceleration 

possible  without  a  skid  occurring  is 

2  2 

M-  g  =  -  0.75  x  32  ft/sec.   =  -  24  ft/sec.  . 
s 

2 

The  leveling  off  of  the  u>  plot  at  -  94  rad/sec.   occurs  at 

the  time  when  the  rotational  velocity  uj  crosses  the  60  rad/sec. 

mark  which  is  the  value  established  for  the  critical  value  of 

rotational  velocity  U)  ,  below  which  u.  becomes  a  constant  =  0.8. 

Since  M*  is  constant  and 

Ff  =  V^mg 

then  the  net  torque  T  is 

T  =  R_Fp,  -  R^^ng  =  a  constant,  and  to  is  thus  a 

constant . 

The  magnitude  of  the  error  rises  to  a  maximum  value  and  then 
decreases  to  zero  when  the  vehicle  has  completely  stopped.   This 
behavior  is  expected  when  it  is  noted  that  the  error  is  merely  the 
difference  between  the  linear  velocity  and  the  rotational  velocity 
of  the  wheel. 

One  other  comment  concerning  the  response  of  the  uncorrected 
system  is  appropriate.   Inspection  of  the  vehicular  velocity  v  plot 
reveals  that  a  steeper  slope  is  present  on  the  plot  prior  to  the 
beginning  of  a  skid.   After  the  skid  begins  the  slope  becomes  less 
steep  and  remains  constant  until  the  vehicle  comes  to  rest.   The 
deceleration  of  the  vehicle  is  thus  greater  before  the  skid  begins 
than  it  is  after  skidding  occurs,  which  is  in  accordance  with  the 
expected  results. 

The  important  fact  to  notice  is  that  during  a  large  part  of 
the  path  the  vehicle  is  in  what  would  be  an  uncontrolled  skid,  as 
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shown  by  the  large  magnitude  of  the  error.   This  uncontrolled  skid 
is  a  matter  of  as  much  concern  as  the  lengthened  time  to  stop  the 
vehicle.   This  lengthened  time  is  caused  by  the  effect  of  the 
smaller  coefficient  of  friction  which  acts  between  the  tire  and 
the  road  surface  when  skidding  is  occurring.   The  optimum  time  of 
stopping,  assuming  braking  takes  place  just  at  the  verge  of  skidding, 

is 

75  ft/sec.      0  1Q 

t  optimum  =  — — - — ; p    =3.13  sec, 

24  ft/sec  ,c 

2 
where  24  ft/sec.   is  the  maximum  possible  deceleration  using  the 

values  assumed  for  this  simulation. 
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5.      THE   DOUBLE   RAMP   CONTROL 

The  various  proposed  controls  described  in  Section  3  were 

implemented  in  an  attempt  to  obtain  a  more  nearly  optimal  braking 

situation.   The  first  of  these  corrective  forces  was  the  double 

ramp  control  shown  in  Fig.  8.   A  total  of  21  separate  runs  were 

made  with  this  form  of  control.   K,  was  initially  held  constant 

at  40000  lbs/sec,  and  K„  was  varied  from  5000  lbs/sec.  up  to 

60000  lbs/sec.   Then  K„  was  held  constant  at  60000  lbs/sec.  and 

K,  varied  from  5000  lbs/sec.  to  50000  lbs/sec.   Finally,  two 

runs  were  made  with  K„  and  K.  set  at  80000  lbs /sec.  and  at  100000 

3      4 

lbs/sec.   Representative  samples  of  the  behavior  to  be  expected 

from  this  form  of  control  are  shown  in  Fig,  15  where  K  =  60000 

lbs/sec.  and  K,  =  35000  lbs/sec,  and  in  Fig.  16  where  K  = 

100000  lbs/sec.  and  K.  =  100000  lbs/sec. 

4 

An  examination  of  Fig.  15  indicates  that  for  the  specific  values 
of  K  and  K,  chosen  (i.e.,  60000  lbs/sec.  and  35000  lbs/sec.  respec- 
tively) the  stop  time  has  been  increased  to  5.47  sec,  whereas  the 
stop  time  under  controlled  conditions  was  only  4.64  sec.   However, 
the  total  time  during  which  a  skid  is  occurring  has  been  reduced  to 
2.08  sec.  and  the  magnitude  of  the  skid  has  been  reduced  to  approx- 
imately 2  ft/sec.  during  the  entire  braking  cycle.   These  figures 
are  contrasted  to  a  total  skid  time  of  4.47  sec.  and  a  maximum  magni- 
tude of  skid  of  68  ft/sec.  for  the  uncorrected  situation. 

When  the  brakes  are  initially  applied  at  time  =  0  sec.  there  is 
a  time  lapse  of  0.16  sec.  before  a  skid  occurs.   This  is  exactly  the 
same  situation  that  occurs  in  the  uncorrected  braking  and  is  accounted 
for  by  the  same  reasoning.  After  0.16  sec,  a  skid  commences.   This 
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is  indicated  in  Fig.  15  by  the  positive  going  error  signal  and  also 

2 
by  the  abrupt  drop  of  the  (U  plot  from  -  24  rad/sec.   to  approximately 

2 
-  78  rad/sec.  .   Recall  that  when  skidding  occurs  the  net  torque 

acting  to  slow  the  wheel  is  suddenly  increased,  thus  the  wheel 

starts  to  slow  down  rapidly,  as  seen  from  inspection  of  the  ou 

versus  time  plot  at  the  time  of  the  initial  occurrence  of  the  first 

skid.   However,  the  corrective  force  immediately  begins  acting  to 

decrease  the  net  force  F  (=  F.   -  F  )  acting  on  the  brake  shoe, 

in    c  ' 

which  consequently  decreases  the  net  torque  T  (=  R-Fp,  -  R  jj,  rng)  , 

and  the  wheel  begins  to  speed  up  to  "catch  up"  with  the  velocity  of 

the  vehicle.   This  accounts  for  the  sudden  positive  slope  of  the  CD 

plot  almost  immediately  after  an  error  signal  occurs,  and  also  for 

the  upswing  on  the  u>   versus  time  plot  after  the  initial  drop  of  u> 

with  each  occurrence  of  a  skid.   The  corrective  force  increases  until 

it  equals  F.   and/or  e  =  0  once  more,  at  which  time  it  begins  decreasing 

towards  zero  until  a  skid  occurs  again. 

An  expanded  plot  of  uo  and  v  versus  time  for  two  skid  cycles  is 

shown  in  Fig.  17.   To  obtain  this  plot  the  equations  of  motion  of  the 

system  were  solved  numerically  by  slide  rule  with  the  value  of  K_  = 

5000  lbs/sec.  and  K,  =  8000  lbs/sec.  Also,  the  simplifying  assumption 

was  made  that  the  coefficient  of  friction  between  the  brake  drum  and 

brake  shoe  was  constant  at  0.8.   While  these  values  for  K«  and  K. 

3  4 

and  the  simplifying  assumption  for  M-  were  not  used  in  the  analog 
computer  runs,  the  shape  of  the  curves  is  similar  to  the  shape  of 
those  obtained  on  the  analog  computer  simulation.   A  brief  explanation 
of  Fig.  17  will  facilitate  understanding  the  curves  shown  in  Fig.  15 
and  Fig.  16. 


38 


•H 
O 

e 

H 
Z 

« 
o 

J3 
• 
► 

1 
I? 

•H 
O 

o 


C 


e 

t 
i 
1 
I 


c 

o 

O 

I 


1 

IS 

u 

t 


-<r 


(oss/q.^    fiqjOO|3/\ 


39 


Referring  to  Fig.  17,  note  that  a  skid  begins  at  time  =  0 

sec.   This  is  due  to  the  simplifying  assumption  of  a  constant 

|J.  =  0.8  over  the  entire  velocity  range  of  interest,  and  in  no  way 

affects  the  qualitative  aspects  of  the  curve  shown.   The  value  of 

uo  (or  R  0)  as  plotted  in  Fig.  17)  decreases  initially  and  then  rises 

to  intersect  the  velocity  V  plot  (which  is  decreasing  linearly  as 

the  vehicle  skids).  For  time  less  than  the  point  marked  (T)  ,  the 

slope  of  the  uu  curve  would  be  positive,  as  can  be  seen  by  mental 

differentiation  of  the  R..UU  curve.   This  is  in  agreement  with  the  tu 

curve  of  Fig.  15  and  Fig.  16.   At  the  point  marked  (l) ,  however, 

the  R-UU  curve  now  begins  increasing  linearly,  and  the  U)  curve  would 

have  a  constant  value,  i.e.,  would  level  off  at  some  positive  value 

corresponding  to  the  slope  of  the  segment  from  (T)  to  (2) .   This 

leveling  off  of  the  w   plot  is  observed  also  in  the  analog  computer 

results  of  Fig.  15  and  Fig.  16.   The  time  period  when  the  leveling 

off  occurs  corresponds  to  the  region  between  (Y)  and  (2J    of  the 

inset  plot  of  F  ,  i.e.,  where  F  has  attained  its  maximum  value 
c  c 

and  the  net  torque  is  now  a  constant.  Since  the  net  torque  is 
constant,  the  value  of  0)  is  also  a  constant,  until  the  skid  is 
completely  eliminated. 

At  point  (2)    the  skid  is  eliminated  and  F  begins  dropping  to 
permit  reapplication  of  the  full  brake  force  F.  .   However,  the  net 
force  F  increases  gradually  from  zero  with  a  slope  K, ,  and  for  a 
significant  time  after  the  elimination  of  the  skid  the  velocity 
is  decreasing  at  an  unacceptably  slow  rate.   This  is  seen  by  an 
inspection  of  the  region  of  the  curve  of  Fig.  17  from  point  (2J  to 
(4)  .  At  point  (3)    the  deceleration  has  finally  increased  to  that 
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2 
value  of  deceleration  experienced  in  a  skid,  i.e.,  -  16  ft/sec.  , 

and  at  point  (4)  it  is  sufficiently  great  to  cause  a  skid  to  occur 

once  more.   The  process  then  repeats  itself. 

The  rather  detailed  qualitative  analysis  undertaken  above  serves 

to  point  up  the  major  defects  inherent  in  the  double  ramp  form  of 

control.   If  the  parameters  K„  and  K,  are  chosen  so  as  to  maximize 

the  time  when  a  skid  is  not  occurring,  e.g.,  to  maximize  the  time 

from  point  (2)    to  point  (4)  of  Fig.  17,  then  the  time  from  point  {2) 

to  (3)  will  also  be  maximized.   From  (2)    to  (3)  the  velocity  is 

decreasing  at  less  than  the  optimum  rate,  and  even  at  less  than 

the  deceleration  experienced  when  in  a  skid.   Thus  it  would  seem 

that  by  minimizing  the  skid  time  with  the  double  ramp  control,  we 

are  inadvertently  maximizing  the  actual  time  to  stop.   This  conclusion 

is  borne  out  by  the  results  of  Fig.  18  and  Fig.  19,  where  the  skid 

times  and  stop  times  for  the  various  values  of  K„  and  K.  simulated 

3      4 

on  the  analog  computer  are  plotted. 

Figure  20  is  a  plot  of  the  magnitude  of  the  skid  observed  versus 
K_,  with  K  fixed  at  40000  lb/sec.   The  magnitude  of  the  skid  after 
it  has  obtained  a  constant  value  is  the  value  that  is  plotted  here. 
It  is  seen  that  the  magnitude  of  the  skid  is  very  dependent  upon  the 
value  of  K„,  which  is  an  expected  result  since  the  larger  the  value 
of  K„,  the  more  rapidly  is  the  skid  corrected.   This  also  accounts 
for  the  decrease  in  the  skid  time  with  increasing  K„  observed  in  Fig. 
18.   Experimental  results  indicated  that  the  magnitude  of  the  skid  is 
unaffected  by  the  value  of  K. ,  which  is  natural  since  the  downward 
going  portion  of  the  double  ramp  control  is  only  acting  when  no  skid 
is  present.   The  stop  time  is  increased,  however,  with  increasing  K„ 
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or  K,  ,  since,  as  mentioned  previously,  the  greater  the  proportion 
of  the  total  stop  time  that  is  free  from  a  skid,  then  the  greater 
will  be  the  actual  time  when  the  deceleration  is  less  than  that 
experienced  when  a  skid  is  present. 

Even  though  the  double  ramp  control  has  serious  defects,  it 
nevertheless  represents  a  decided  improvement  over  the  completely 
uncontrolled  braking  represented  in  Fig.  14.   Since  under  many  road 
conditions  it  is  the  loss  of  control  which  occurs  when  a  skid  is 
present  that  is  actually  more  dangerous  than  the  increase  in  the 
stop  time,  then  the  double  ramp  control  would  have  certain  merits. 
The  double  ramp  control  is  capable  of  minimizing  the  skid  time  and 
the  magnitude  of  the  skid,  at  a  cost  of  slightly  increased  stopping 
distance. 

If  it  is  assumed  that  the  stop  time,  the  skid  time,  and  the 
degree  of  skid  experienced  are  of  equal  importance  in  judging  the 
overall  goodness  of  a  particular  type  of  control,  then  a  performance 
factor  (P.F.)  can  be  defined  for  the  uncorrected  case,  the  double 
ramp  control,  and  for  subsequent  controls  to  be  tried.   This  per- 
formance factor  is  defined  as  the  reciprocal  of  the  product  of  the 
actual  stop  time  minus  the  optimum  stop  time,  the  skid  time,  and 
the  degree  of  skid  experienced,,   The  degree  of  skid  will  be  measured 
after  the  velocity  has  decreased  below  60  ft/sec,  since  that  is  the 
velocity  below  which  the  value  for  p.  (between  the  brake  drum  and 
brake  shoe)  is  a  constant.   The  optimum  stop  time  is  3.13  sec. 

This  definition  of  a  performance  factor  must  not  be  taken  to 
be  the  absolute  criterion  for  judging  a  control  and  can  only  be 
used  as  a  rough  guide  to  determine  the  relative  merits  of  different 
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types  of  controls.   It  cannot,  for  example,  be  used  to  determine 
the  proper  setting  of  control  parameters  to  give  the  ideal  stopping 
conditions,  since  the  optimum  settings  for  one  set  of  road  conditions 
might  very  possibly  be  different  from  the  optimum  settings  for  another 
set  of  road  conditions.  Also,  no  account  is  taken,  in  the  definition 
of  a  performance  factor,  of  the  fact  that  a  slight  degree  of  skidding 
might  be  much  less  significant  than  a  larger  degree  of  skidding,  or 
that  an  increase  in  the  stop  time  might  be  permissible  under  some 
circumstances  but  entirely  unacceptable  under  others.   Nevertheless, 
the  definition  given  above  is  convenient  and  useful,  within  limits. 

For  the  uncorrected  situation  of  Fig.  14,  a  performance  factor 
can  be  computed  in  order  to  provide  a  basis  for  comparison.   The 
integral  of  the  error  curve,  rather  than  the  product  of  the  error 
times  the  skid  time,  is  used,  multiplied  by  the  total  stop  time  minus 
the  optimum  stop  time.   This  gives  a  performance  factor  of  0.00214 
for  the  completely  uncontrolled  situation. 

For  the  double  ramp  control,  the  values  of  the  performance 
factor  range  from  0.0154  to  0.298,  and  the  variation  of  P.F.  with 
a  change  in  K„  is  shown  in  Fig.  21.   This  figure  is  presented  only 
to  indicate  the  trend  of  P.F.  for  the  double  ramp  control  and  not  as 
an  aid  to  selecting  the  "best"  value  for  K~ .   However,  based  on  the 
indication  given  by  the  performance  factors  for  the  double  ramp 
control,  the  double  ramp  control  provides  a  definite  overall  perform- 
ance of  the  system,  when  compared  to  the  uncontrolled  situation. 


46 


jo- 


t 


.08  - 


a 
o 

u 

<  .06 
L. 


cj 
< 

5  -of  ■ 

q: 

o 

u. 

or 

u 

a 

.o£- 


.01 


i      i     i      i  i  i  i  i 

\000O       S.OQOO        3O0O0      LjOOOO       50OO0        £00OO  L>       ('Ita/ger)- 

Fig.  21.  Perforaanco  factor  rorsus  K~  for  the  doublo 

ramp  control. 


47 


6.   THE  BANG- BANG  CONTROLS. 

The  bang-bang  control  with  a  dead  zone  was  simulated  on  the 
analog  computer.   The  dead  zone  was  arbitrarily  set  at  7.5  ft/sec., 
which  is  rather  high  for  a  physical  system,  but  which  nevertheless 
serves  to  illustrate  the  behavior  of  the  system  under  this  form 
of  control.   The  response  obtained  is  shown  in  Fig.  22. 

The  stop  time  is  seen  to  be  4.50  sec,  which  is  slightly  less 

than  that  obtained  for  the  completely  uncontrolled  situation.   This 

is  a  predictable  result,  since  for  almost  the  entire  cycle  the 

vehicle  is  in  a  "controlled"  skid,  where  the  term  "controlled" 

implies  that  the  value  of  the  skid  is  never  permitted  to  exceed  a 

certain  value.   The  reason  for  the  slightly  lower  stop  time  than 

that  observed  for  the  uncontrolled  situation  is  that  for  a  small, 

non-zero  time  interval  in  the  region  where  F  drops  to  zero  and  F 

rises  to  F.  ,  the  value  of  the  skid  is  essentially  zero  (or  at  any 
in 

rate  less  than  the  threshold  established  for  detection  of  an  error 
on  the  basic  system  of  Fig.  7).   While  the  skid  is  this  small,  the 
coefficient  of  friction  controlling  the  deceleration  of  the  vehicle 
is  \i,      (which  is  greater  than  p.  ) ,  and  the  deceleration  is  thus 
greater  for  this  small  time  interval.   Since  this  phenomenon  occurs 
every  time  the  skid  is  corrected,  the  cumulative  effect  over  the 
entire  stopping  distance  is  to  decrease  the  stop  time  slightly. 
The  major  disadvantage  to  this  form  of  control  is  the  fact 
that  a  skid,  even  though  of  a  small  value,  is  almost  always  present 
until  the  vehicle  comes  to  a  complete  stop.   This  also  implies  that 
the  shortest  possible  stop  time  is  nearly  4.64  sec,  i.e.,  the  same 
time  required  to  stop  when  no  corrective  force  is  present.   However, 
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the  fact  that  the  skid  is  controlled,  and  can  theoretically  be 
reduced  to  an  arbitrarily  small  positive  value  (by  reducing  the 
dead  zone)  indicates  that  even  this  form  of  control  is  superior  to 
the  uncontrolled  situation.  An  additional  advantage  is  that  as 
the  degree  of  skid  is  reduced  toward  0,  i.e.,  as  the  relative 
velocity  between  the  tire  and  road  surface  is  reduced,  the  coef- 
ficient of  friction  becomes  somewhat  greater  than  that  present 
when  a  more  significant  relative  velocity  is  attained.   This  means 
that  a  greater  value  of  deceleration  is  possible  when  the  skid 
becomes  extremely  small.   This  feature  of  the  coefficient  of 
friction  was  not  incorporated  in  these  analog  computer  simulations 
due  to  the  greatly  increased  complexity  of  the  simulation  that  would 
be  necessitated. 

Even  though  this  control  could  theoretically  provide  good  perform- 
ance if  the  dead  zone  were  reduced  sufficiently,  there  is  a  serious 
physical  limitation  that  prevents  this  from  ever  really  being  possible. 
A  little  reflection  will  reveal  that  the  lower  the  dead  zone  setting, 
the  more  rapidly  will  the  corrective  force  be  required  to  go  on  and 
off.  As  the  dead  zone  approaches  zero,  the  shape  of  the  corrective 
force  will  appear  as  a  high  speed  oscillation,  and  the  practical  lower 
limit  for  the  dead  zone  will  be  determined  by  both  the  speed  of  response 
of  the  hydraulic  or  pneumatic  lines  and  the  physical  strength  of  the 
brake  system. 

Thus,  from  a  practical  point  of  view,  this  system  appears  to  offer 
the  advantage  of  a  controlled  skid,  plus  a  shorter  stop  time  than  that 
provided  by  the  double  ramp  control.   The  disadvantage  is  that  a  skid 
is  almost  always  present,  and  the  level  to  which  it  could  be  reduced 
would  be  severely  limited  by  the  physical  response  of  the  brake  lines. 
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Neither  control  tried  thus  far  appears  to  be  superior  to  the 
other.   This  is  verified  somewhat  by  determining  the  performance 
factor  for  the  bang-bang  control  with  a  dead  zone.   The  performance 
factor  (computed  by  taking  the  reciprocal  of  the  product  of  the 
integral  of  the  error  and  the  difference  between  the  stop  time  and 
the  optimum  time  to  stop)  is  determined  to  be  0.046.   This  is 
neither  significantly  better  nor  worse  than  any  of  the  performance 
factors  obtained  for  the  double  ramp  control. 

The  bang-bang  control  with  a  time  delay  was  then  simulated. 
The  results  for  this  simulation  are  presented  in  Fig.  23.  A  time 
delay  of  0.1  sec.  (problem  time)  was  chosen  for  this  simulation. 
The  stop  time  is  seen  to  be  5.93  sec,  which  is  considerably  worse 
than  the  stop  time  for  the  previously  described  bang-bang  control. 
However,  in  this  case  the  skid  is  eliminated  for  a  portion  of  the 
total  stop  cycle,  and  the  skid  time  is  3.92  sec„   Superficially  this 
appears  to  be  an  improvement  over  the  previous  bang-bang  control, 
but  closer  examination  reveals  that  very  little,  if  any,  betterment 
of  the  performance  is  gained  by  the  incorporation  of  a  time  delay. 

Reference  to  Fig*  24,  which  is  an  expansion  of  a  portion  of 
Fig.  23,  and  was  obtained  by  manual  solution  of  simplified  equations 
of  motion,  reveals  that  for  a  significant  period  of  time  during  each 
error  correcting  cycle  the  velocity  remains  a  constant.   For  example, 
from  point  f2~)  to  (J)  ,  where  no  skid  is  present,  the  velocity  is 
also  a  constant.   The  time  from  (2)     to  (3)  is  0.1  sec,  and  corresponds 
to  the  time  delay  built  into  the  corrective  force.   During  this  time, 
the  net  force  F  on  the  brake  shoe  is  zero,  and  thus  the  velocity 
remains  constant.   If  the  time  delay  is  decreased,  this  control 
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approaches  the  performance  of  the  bang-bang  control  with  a  zero 
dead  zone. 

The  bang-bang  control  with  a  time  delay  thus  appears  superior 
to  a  completely  uncorrected  braking  situation,  but  appears  some- 
what inferior  to  either  of  the  two  controls  previously  described. 
A  performance  factor  for  this  control  (for  the  specific  set  of 
parameters  used  for  this  simulation)  is  0.025,  which  is  sufficiently 
lower  than  the  trend  of  the  values  obtained  for  the  previous  two 
controls  to  warrant  the  above  conclusion. 
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7.   THE  TRIPLE  RAMP  CONTROL. 

The  triple  ramp  control  of  Fig„  12  was  the  final  type  of 
corrective  scheme  to  be  implemented .   A  total  of  100  separate 
computer  runs  were  made,  using  a  variety  of  values  for  the 
parameters  K  ,  K  ,  K  and  C.R.  (correction  ratio).   Representative 
samples  of  the  behavior  to  be  expected  from  this  type  of  corrective 
scheme  are  shown  in  Fig.  25  and  Fig.  26. 

The  most  striking  features  about  both  Fig.  25  and  Fig„  26 

are  the  low  stop  time,  skid  time,  and  magnitude  of  error  present 

for  this  form  of  control.   For  example,  Fig.  25,  with  K„  =  60000 

lbs/sec,  K.  =  40000  lbs/sec,  K  =  1000  lbs/sec.  and  C.R.  =  2.05 
4  c 

has  a  stop  time  of  4.22  sec.,  skid  time  of  0.98  sec.  and  roughly 
1.9  ft/sec.  as  the  degree  of  skid.   This  gives  a  performance  factor 
of  0.493  which  is  higher  than  the  maximum  value  attained  for  any 
control  tried  previously.   This  value  of  performance  factor  is  of 
course  affected  by  the  values  of  the  various  parameters  K„ ,  K,  ,  K  , 
and  C.R.,  and  does  not  represent  the  "best"  value  of  P.F.  that  was 
obtained  on  the  trial  runs  for  the  triple  ramp  control.   Indeed, 
the  term  "best"  is  misleading  here,  since  a  setting  of  parameters 
resulting  in  the  highest  possible  P.F.  would  sacrifice  an  adapta- 
bility to  changing  road  conditions.   This  feature  will  be  discussed 
more  thoroughly  later  in  this  paper. 

Nevertheless,  it  can  be  concluded  from  the  above  factors  and 
an  inspection  of  Fig.  25  and  Fig.  26  (along  with  the  Author's 
analysis  of  the  other  98  analog  computer  runs)  that  the  triple  ramp 
control  provides  more  nearly  optimal  braking  conditions  than  any 
previously  described  control. 
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A  closer  examination  of  the  actual  behavior  of  this  triple 

ramp  control  is  appropriate  at  this  point.   The  results  of  a 

manual  solution  of  the  equations  of  motion  (with  the  simplifying 

assumptions  mentioned  previously  for  similar  hand  calculations)  are 

presented  graphically  in  Fig.  27.   The  parameter  values  assumed 

were  K0  =  50000  lbs/sec,  K,  =  80000  lbs/sec,  K  =  1000  lbs/sec, 
3  4  c 

and  C.R.  =  1.5.  Also,  the  coefficient  of  friction  ijl  between  the 
brake  drum  and  brake  shoe  was  assumed  to  be  constant  at  0.8  for 
the  hand  calculation  (this  was  the  only  simplifying  assumption). 
All  other  conditions  were  assumed  the  same  as  previously  listed 
for  other  simulations. 

An  examination  of  Fig.  27  indicates  that  from  time  =  0  sec. 
to  the  point  marked  (5)  ,  the  behavior  of  the  control  and  the 
vehicle  response  are  precisely  analogous  to  the  double  ramp  form 
of  control.   This  can  best  be  seen  by  a  comparison  of  Fig.  27 
with  Fig.  17  over  the  region  from  time  =  0  sec.  to  the  point  (5)     . 
The  same  defect  present  in  the  double  ramp  control  is  also  present 
in  the  triple  ramp  control,  i.e.,  from  point  (Z)    to  point  \3) 
the  velocity  is  decreasing  at  a  rate  slower  than  the  rate  of 
deceleration  when  a  skid  is  present.  Also,  from  point  (5y  to  point 
(6)    this  phenomenon  occurs.  However,  at  point  (7)  ,  when  the  down- 
ward slope  of  the  corrective  force  F  abruptly  changes  from  -  K, 
to  -  K  the  situation  changes  completely.  While  the  slope  of 
F  =  -  K  ,  the  rate  of  decrease  of  the  velocity  not  only  is  signif- 
icantly greater  than  that  possible  during  a  skid,  but  also  no  error 
is  present  until  the  value  of  the  deceleration  again  exceeds  -  24.0 
ft/sec.   at  point  (s)  .  At  this  time  the  process  repeats  itself. 
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It  is  during  this  last  phase  of  the  skid  correcting  cycle 

(from  point  (j)    to  point  (§)  in  Fig.  27)  that  the  inherent 

superiority  of  the  triple  ramp  control  is  manifested.   During 

this  portion  of  the  cycle  the  deceleration  is  gradually  increasing 

as  F   slowly  decreases.   When  the  deceleration  is  too  great  for 
c 

the  prevailing  road  conditions  a  skid  once  again  occurs,  but  prior 
to  this  occurrence  of  the  skid  the  deceleration  is  more  and  more 
closely  approaching  the  optimum  deceleration.  Additionally,  if 
the  road  conditions  are  changing,  for  example,  if  the  vehicle 
went  from  a  patch  of  slick  pavement  to  dry  pavement,  and  the  coef- 
ficient of  friction  between  the  tire  and  road  surface  were  increased, 
then  this  form  of  control  permits  the  deceleration  to  increase  to 
the  new  optimum  value.   In  other  words,  the  triple  ramp  control 
is  continually  sampling  the  road  conditions  and  permitting  the 
braking  force  to  adapt  to  changing  road  conditions. 

The  desirability  of  this  adaptability  feature  is  evident.   It 
is  also  apparent  that  for  the  invariant  road  conditions  simulated 
for  this  problem,  the  most  nearly  ideal  braking  situation  would 
have  been  achieved  by  setting  C.R.  to  a  value  just  slightly  greater 
than  one,  and  causing  K  to  equal  zero.   This  would  have  permitted 
only  two  skids  of  small  magnitude  to  occur,  and  for  the  remainder 
of  the  path  the  vehicle  would  have  slowed  down  at  a  rate  very  nearly 
equal  to  the  optimum  deceleration.   This  would  provide  a  very  large 
value  for  the  performance  factor,  which  would  of  course  be  a  misleading 
indication  in  this  case,  since  the  control  would  not  respond  to  improved 
road  conditions.   It  would,  however,  still  respond  to  worsening  road 
conditions,  so  there  is  even  some  merit  in  this  proposal. 
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Referring  back  to  Fig.  25,  we  note  that  the  typical  triple 
ramp  shape  of  F  does  not  become  evident  until  after  the  occurrence 
of  the  fourth  error  signal,  whereas  in  Fig.  26  the  triple  ramp 
shape  appears  as  expected  after  the  occurrence  of  the  second  error 
signal.   Recall  that  the  form  of  the  coefficient  of  friction  |i 
between  the  brake  drum  and  brake  shoe  is  given  by 

[i,   =    .8  -  .01  (u)  -  60)  for  ud  >  60  rad/sec. 

Reference  to  the  uo  versus  time  plot  of  Fig.  25  reveals  that  both 

the  first  and  second  errors  occur  while  the  rotational  velocity  is 

above  60  rad/sec.   Thus  in  this  region  the  value  for  fJb  is  constantly 

changing  in  accordance  with  the  above  relation.   Now,  when  the  first 

error  signal  (or  skid)  is  sensed  the  corrective  force  F  immediately 

increases  until  the  skid  is  eliminated,  and  then  decreases  until 

the  skid  occurs  once  more;  this  value  of  F  at  which  the  skid  occurs 

'  c 

is  stored  by  a  storage  device.   Upon  the  occurrence  of  the  second 

skid,  F   increases  again  until  the  skid  is  once  more  eliminated, 
c 

Ideally  the  value  of  F  should  drop  to  some  multiple  of  the  stored 
previous  value  of  F   (the  multiple  being  established  by  C.R.).   How- 
ever, the  velocity  of  the  vehicle  is  decreasing  continually  during 
this  process  and  [i.   is  thus  increasing.   This  means  that  after  the 
occurrence  of  the  second  skid,  the  net  torque  on  the  wheel  for  a 
given  value  of  F   (e.g.,  for  a  value  equal  to  C.R.  times  the  previously 
stored  value  of  F  )  is  now  greater  than  it  would  be  if  the  value  for  (J. 
were  not  increasing.   Thus  a  third  skid  occurs  before  F  can  drop  to 
the  value  at  which  the  break  point  occurs.   However,  after  the  occur- 
rence of  the  fourth  skid,  p.  is  constant,  and  thus  the  typical  triple 
ramp  effect  is  observed. 
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In  Fig.  26  the  above  phenomenon  does  not  occur.   Note  that 

in  this  figure  the  value  of  K-  is  only  25000  lbs /sec.  whereas  in 

Fig.  25  K„  had  a  value  of  60000  lbs/sec.   This  means  that  the 

error  will  be  corrected  more  slowly,  and  thus  the  velocity  will 

have  decreased  to  a  level  where  p,  is  a  constant  by  the  time  the 

second  error  has  been  eliminated.   Thus  the  coefficient  of  friction 

will  not  have  changed  sufficiently  between  the  time  when  the  value 

of  F  was  stored  and  the  time  when  F  reaches  its  break  point  to 
c  c  r 

cause  another  error  signal  to  occur.   Note,  however,  that  if  C.R. 

had  been  set  sufficiently  low  (e.g.,  slightly  above  1.0)  then  the 

same  phenomenon  that  occurred  in  Fig.  25  would  also  have  occurred 

in  Fig.  26.   Conversely,  had  C.R.  been  set  sufficiently  high,  this 

phenomenon  would  not  have  occurred  in  Fig.  25. 

Figure  28  presents  a  plot  of  stop  time  versus  K  ,  with  C.R. 

as  a  parameter.   The  data  points  on  this  and  subsequent  plots 

were  obtained  from  the  graphical  output  resulting  from  the  various 

analog  computer  runs  made  while  simulating  the  triple  ramp  control. 

Consequently,  the  accuracy  of  the  individual  data  points  is  only 

of  the  order  of  a  few  percent.   However,  this  degree  of  accuracy  is 

sufficient  to  demonstrate  the  overall  performance  characteristics 

of  the  triple  ramp  control  quite  satisfactorily.   From  Fig.  28  we 

can  see  that  as  K  increases  the  stop  time  increases  also.   However, 
c 

the  rate  of  increase,  i.e.,  the  slope  of  the  curve,  is  observed  to 
be  a  variable,  even  though  the  overall  trend  is  upward.   This  over- 
all upward  trend  is  expected  in  view  of  previous  comments  concerning 
the  adaptability  feature  inherent  in  the  triple  ramp  control.  A 
decrease  in  stop  time  must  be  achieved  at  the  cost  of  a  decrease  in 
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the  sensitivity  of  the  control  to  improving  road  conditions,  all 
other  parameters  being  assumed  constant.   Thus,  as  K  increases, 
the  vehicle  goes  into  a  skid  mode  more  frequently,  and  the  overall 
trend  of  stop  time  with  K  would  be  an  upward  one. 

However,  to  satisfactorily  account  for  the  leveling  off  of 
the  stop  time  versus  K  chart  at  localized  regions,  we  must  investi- 
gate the  detailed  behavior  more  closely.   Figure  29  is  a  plot  of 
both  stop  time  and  skid  time  versus  K  for  one  specific  set  of  para- 
meters.  The  numbers  adjacent  to  the  skid  time  curve  represent  the 
total  number  of  skids  present  in  the  particular  graphical  output 
from  which  each  specific  data  point  was  obtained.   For  example,  8 
skids  were  present  in  the  analog  computer  output  obtained  with  a 
setting  of  K  =  60000  lbs/sec,  K  =  25000  lbs/sec,  C.R.  =  2.05, 

and  K  =  1500  lbs/sec.   For  K  =  2500  lbs/sec.  and  K  =  3000  lbs/sec, 
c  c  c 

10  skids  were  present  and  the  skid  time  was  thus  the  same  for  each 

of  these  values.   However,  Fig.  29  also  reveals  that  the  stop  times 

for  each  of  these  values  of  K  were  the  same. 

c 

In  a  situation  where  the  same  number  of  skids  is  observed  from 

one  value  of  K  to  another,  then  the  skid  times  will  obviously  be 
c 

identical,  assuming  that  K„  (which  directly  determines  the  magnitude 
of  each  skid  and  the  length  of  time  it  is  present)  is  the  same  for 
each  trial.   The  only  factor  contributing  to  a  difference  in  over- 
all stop  time  is  thus  the  difference  in  time  between  adjacent  pairs 
of  skids  (e.g.,  the  time  from  point  (j5)  to  point  (&)     in  Fig. m  27). 
For  a  larger  value  of  K  ,  the  time  between  pairs  of  skids  is  decreased. 
Now,  for  example,  when  10  skids  occur  for  a  value  of  K  .=  3000  lbs/sec. 
the  time  between  adjacent  pairs  of  skids  is  shorter  than  the  time 
between  the  adjacent  pairs  of  10  total  skids  obtained  for  K  =  2500 
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lbs/sec.   Thus  the  final  pair  of  skids  occurs  earlier  when  K  = 

c 

3000  lbs/sec,  but  due  to  a  larger  value  of  K  a  larger  value  of 

deceleration  is  attained  after  the  tenth  skid  and  prior  to  the 

vehicle  coming  to  a  complete  stop.   There  is  also  a  longer  period 

of  time  after  the  tenth  skid  occurs,  during  which  the  slope  of 

F  =  -  K  ,  for  the  case  where  K  =  3000  lbs/sec.   These  two 
c      c  c 

factors  serve  to  account  for  the  leveling  off  of  the  stop  time 
versus  K  curves  observed  for  Fig.  28  and  Fig.  29.   It  is  also 
possible  that  a  slight  local  decrease  in  total  stop  time  might 
occur  for  an  increase  in  K  under  the  proper  combination  of 
control  parameters  and  road  conditions.   This  is  observed  in  Fig. 
28  in  one  region  of  the  curve  obtained  for  C.R.  =  2.74. 

Despite  the  above,  observations,  the  conclusion  to  be  drawn 
from  Fig.  28  and  Fig.  29  is  that  the  overall  effect  of  an  increase 
in  K  is  to  increase  the  stop  time  of  the  vehicle,  and  to  degrade 
the  overall  performance  of  the  corrective  system.   A  compromise 
must  thus  be  reached  between  adaptability  to  improving  road 
conditions  and  decreasing  the  total  stop  time.   Designing  for  a 
worst  case  situation,  a  satisfactory  value  for  K  would  thus  be 
K  =  500  lbs/sec,  as  this  would  permit  reasonably  rapid  response 
to  improving  road  conditions  while  still  providing  very  good  braking 
performance  under  nearly  static  road  conditions. 

The  next  parameter  to  be  investigated  is  the  correction  ratio 
(C.R.).   Figure  30  is  a  plot  of  the  stop  time  versus  the  correction 
ratio,  with  K  as  a  parameter.   The  overall  trend  as  determined  from 
Fig.  30  is  that  the  stop  time  tends  to  decrease  with  an  increase  in 
the  value  of  C.R.   Each  curve  in  the  family  of  Fig.  30  demonstrates 
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notable  local  exceptions  to  this  trend,  however,  and,  as  will 

be  pointed  out,  excessively  large  values  of  C.R.  will  materially 

lengthen  the.  stop  time. 

The  overall  downward  trend  of  the  stop  time  versus  C.R.  curve 

is  best  understood  by  referring  to  the  inset  of  Fig.  27,  where  a 

pictorial  description  of  the  triple  ramp  control  is  presented. 

As  the  value  of  C.R.  is  raised  this  raises  the  break  point  at 

which  the  transition  of  the  slope  of  F  from  -  K.  to  -  K  occurs. 

c         4       c 

This  has  the  effect  of  lengthening  the  time  between  the  occurrence 

of  adjacent  skid  pairs,  since  the  time  that  F  will  take  to  drop 

to  the  level  where  a  skid  will  occur  once  more  increases  in  direct 

proportion  to  the  height  of  the  break  point.   It  is  thus  seen  that 

increasing  C.R.  has  the  same  net  effect  as  decreasing  the  value 

of  K  ,  with  certain  minor  restrictions, 
c 

It  would  be  expected  that,  as  C.R.  is  increased  to  obtain  a 

set  of  braking  histories,  those  output  traces  which  have  the  same 

number  of  skids  present  for  differing  values  of  C.R.  would  behave 

in  a  manner  similar  to  that  observed  when  K  was  the  variable 

c 

quantity  (e.g.,  in  Fig.  28  and  Fig.  29).   Referring  in  Fig.  30  to 

the  curve,  obtained  for  K  =  500  lbs/sec.,  we  see  that  the  stop  time 

c 

increases  with  increasing  C.R.  when  the  total  number  of  skids  present 

in  each  braking  history  remains  a  constant  four  (this  is  also  observed 

for  K  =  1000  lbs/sec.,  for  six  skids  observed).   The  reason  for 
c 

this  is  very  simple.  As  C.R.  is  raised  this  "pushes"  the  adjacent 
pairs  of  skids  farther  apart  and  thus  for  a  larger  portion  of  the 
total  stop  time  the  vehicle  is  decelerating  at  a  lower  rate  than  if 
C.R.  were  set  lower.   Thus  if  the  total  of  skids  still  remained  the 
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same,  the  control  which  provided  the  greater  deceleration  over 

the  entire  skid  free  path  would  also  provide  a  shorter  stop  time. 

This  is  the  control  with  the  lower  C.R. 

Note  that  C.R.  is  also  closely  linked  to  the  adaptability  of 

the  system  to  improving  road  conditions,  since  for  a  large  value 

of  C.R.  the  time  required  for  F  to  drop  (at  a  rate  of  -  K  )  to 

c  c 

a  level  where  the  road  is  "tested"  for  better  conditions  is 
lengthened.   Additionally,  C.R.  cannot  be  increased  to  a  very 
large  value  without  greatly  extending  the  stop  time.   The  minimum 
number  of  skids  possible  (starting  with  any  significant  initial 
velocity)  is  two.   Once  C.R.  has  been  increased  to  a  level  where 
only  two  skids  are  present,  then  any  further  increases  in  C.R.  will 
cause  a  rapid  increase  in  the  stop  time  for  the  reasons  listed  in 
the  preceding  paragraph.   However,  long  before  the  skids  were 
reduced  to  two  in  number,  the  stop  time  would  also  have  increased 
due  to  the  low  value  of  deceleration  that  would  be  present  over  an 
excessively  large  portion  of  the  vehicle's  trajectory. 

In  retrospect,  it  appears  that  despite  the  downward  trend  of 
stop  time  with  increasing  C.R.  presented  in  Fig.  30,  any  advantages 
to  be  gained  by  very  large  values  of  C.R.  are  outweighed  by  the 
disadvantages.   On  the  basis  of  information  presented  here,  a  "safe" 
value  for  C.R.  would  appear  to  be  about  2.0. 

Figure  31  and  Fig.  32  present  plots  of  both  stop  time  and  skid 
time  versus  C.R.  for  two  different  sets  of  parameters.   They  are 
included  primarily  for  information,  the  most  significant  factor 
about  either  one  being  that  where  the  curve  representing  skid  time 
becomes  flat  (implying  a  constant  number  of  skids  were  observed  for 
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succeeding  braking  history  plots)  the  curve  for  stop  time  acquires 
a  positive  slope.   This  serves  to  reinforce  the  arguments  previously 
advanced. 

The  behavior  of  the  braking  performance  as  K„  is  varied  is 
also  of  some  interest.   Figure  33  illustrates  the  behavior  of  stop 
time  and  skid  time  as  K  is  increased.   Notice  that  the  tendency 
exhibited  is  for  the  stop  time  to  decrease  with  increasing  K_.   This 
is  the  expected  result,  since  the  larger  the  value  of  K_,  the  more 
rapidly  is  any  error  or  skid  corrected  back  to  zero,  and  thus  a 
larger  portion  of  the  braking  path  will  be  in  a  controlled  mode  of 
skid  free  braking. 

Referring  back  to  Fig.  27  once  more,  we  note  that  the  time 
from  the  commencement  of  braking  action  (t  =  0  )  to  point  (T)  ,  and 
also  from  point  (4)  to  point  £3)  ,  is  dependent  solely  upon  the  value 
of  K_,  whereas  the  time  from  point  (2)  to  point  (_4)  a^d  from  point  (5) 
to  point  (8)  depends  upon  K,  and  K  .   Thus  as  K  is  increased  the 
total  skid  time  and  magnitude  of  each  individual  skid  are  reduced. 
This  has  the  effect  of  causing  each  successive  skid  pair  to  occur 
earlier  in  time,  even  though  the  time  between  skids  is  unchanged. 
As  long  as  the  total  number  of  skids  present  is  unchanged  as  K„  is 
increased,  an  improvement  will  always  be  observed  in  the  stop  time. 
However,  in  the  event  that  one  value  of  K„  causes  a  certain  number 
of  skids,  and  a  higher  value  of  K„  is  just  sufficient  to  increase 
the  number  of  skids,  then  a  slight  local  increase  will  be  observed 
in  the  stop  time.   This  phenomenon  is  seen  to  occur  in  Fig.  33. 
The  explanation  for  this  is  as  follows.  When  K  =  50000  lbs/sec. 
the  time  remaining  between  the  correction  of  the  sixth  skid  and  the 
complete  stopping  of  the  vehicle,  is  insufficient  for  the  deceleration 
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to  have  increased  to  that  level  necessary  to  cause  another  skid. 
However,  for  K„  =  60000  lbs/sec.  the  seventh  skid  occurs  just 

prior  to  the  vehicle  coming  to  a  rest.   During  this  seventh  skid 

2 
the  vehicle  deceleration  is  only  -  16  ft/sec.  ,  rather  than  the 

2 
-  24  ft/sec.   that  was  experienced  just  prior  to  the  skid.   It  is 

this  lower  value  of  deceleration  present  just  prior  to  stopping 

that  causes  the  stop  time  to  increase  as  K~  is  increased  from 

50000  lbs/sec.  to  60000  lbs/sec. 

The  above  phenomenon  is  just  restricted  to  the  transition 
points  where  the  number  of  complete  skids  increases  from  one  value 
to  another.   However,  the  decreased  skid  times  and  magnitude  of 
skids  caused  by  an  increased  value  of  K„  indicates  that  the  greatest 
overall  advantage  is  to  be  gained  by  making  K„  as  large  as  possible, 
regardless  of  any  localized  increases  in  stop  time  caused  by  this 
policy. 

Finally,  the  system  behavior  with  changing  values  of  K,  is 
investigated.   Figure  34  indicates  that  the  stop  time  decreases  as 
the  value  of  K,  is  increased,  but  that  the  skid  time  goes  up  with 
increasing  K, .   Reference  to  Fig.  27  again  facilitates  understanding 
this  behavior.   From  Fig.  27  it  can  be  seen  that  the  time  lapse  from 
point  (2)    to  point  (4/  and  also  from  point  (j)  to  point  (j)     is  depend- 
ent solely  upon  the  value  of  K, .   Large  values  of  K,  will  cause  the 
time  intervals  between  the  points  to  become  small,  while  leaving 
time  intervals  along  other  portions  of  the  trajectory  unchanged. 
Since  it  is  during  the  time  interval  from  point  (5)  to  point  (§)  , 
for  example,  that  the  deceleration  of  the  vehicle  is  even  less  than 
that  present  while  in  a  skid,  any  factor  which  will  shorten  the  time 
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from  (5)  to  (j)    will  also  shorten  the  time  the  vehicle  spends  in 
the  undesirable  region  from  (5)  to  (6)    .   This  will  naturally 
cause  a  shorter  stop  time,  since  the  vehicle  will  be  decelerating 
at  a  high  rate  over  a  larger  proportion  of  its  stopping  path. 

The  increase  of  skid  time  noted  in  Fig.  34  is  once  again 
related  closely  to  the  total  number  of  skids  present.  As  the  time 
from  (T)    to  (5)  or  from  (5)  to  (7)  (referring  once  again  to  Fig.  27) 
is  reduced,  this  tends  to  shorten  the  time  between  the  occurrence 
of  adjacent  skids.   If  the  total  number  of  skids  remains  the  same 
while  K,  is  being  increased,  the  skid  time  naturally  remains  the 
same  (to  fully  appreciate  this,  one  must  recall  that  the  duration 
and  magnitude  of  each  individual  skid  is  controlled  by  K„  only, 
which  is  not  being  varied  here).   However,  as  the  skid  occurrences 
become  closer  together  in  time,  a  point  is  reached  where  the  total 
number  of  skids  encountered  in  the  stopping  history  increases,  and 
thus  the  skid  time  will  also  increase.   This  increase  in  skid  time 
would  not  be  a  serious  consequence,  however,  if  K„  were  sufficiently 
high  to  make  the  duration  and  magnitude  of  the  individual  skids  very 
sma 1 1 . 

Thus  it  might  be  concluded  that  the  larger  the  value  of  K  , 
the  better  the  performance  of  the  entire  system,  since  the  decrease 
in  stop  time  would  more  than  offset  any  increase  in  skid  time,  partic- 
ularly if  K»  were  also  set  very  high.   Caution  must  be  exercised  here, 
however,  because  it  is  the  nature  of  the  triple  ramp  control  that 
the  value  of  F  at  which  every  second  error  occurs  must  be  marked 
and  stored.   This  means  that  if  K,  is  too  large,  then  difficulties 
would  be  encountered  in  marking  the  point  at  which  each  second  error 
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occurs,  and  additionally  in  reversing  the  trend  of  F  (i.e.,  in 
changing  from  a  slope  of  -  K,  to  +  K_).  Obviously  a  slope  of 
-  infinity  for  -  K,  would  be  totally  unsatisfactory.  Thus,  a 
very  large  value  of  K,  would  indeed  provide  the  best  performance, 
but  this  value  must  not  be  so  great  as  to  create  excessive  dif- 
ficulties in  measuring  and  storing  the  point  at  which  each  second 
error  occurs . 

In  summary,  the  following  qualitative  comments  concerning 
the  choice  of  the  values  for  the  various  parameters  of  the  triple 
ramp  control  are  appropriate.  K»  should,  with  no  qualifications, 
be  made  as  large  as  possible.  K,  should  also  be  as  large  as  possible, 
short  of  being  a  negative  step  function,  and  subject  to  the  restric- 
tions mentioned  in  the  preceding  paragraph.  K  should  be  relatively 
low  assuming  that  C.R.  is  not  too  large.  A  value  for  C.R.  of  about 
2.0,  and  for  K  of  500  lbs/sec.  was  determined  to  be  satisfactory 
for  the  particular  system  simulated  here. 
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8.   PROPOSED  IMPLEMENTATION  OF  THE  TRIPLE  RAMP  CONTROL . 

In  view  of  the  obvious  superiority  of  the  triple  ramp  form 
of  control  over  any  of  the  other  controls  examined,  a  proposal 
for  the  physical  implementation  of  this  control  is  advanced  here. 

It  must  be  assumed  that  a  linear  accelerometer  of  high  accuracy 

2  2 

and  linearity  over  the  range  from  0  ft/sec.   to  32  ft/sec.   is 

available o  A  high  gain,  broad  bandwidth  operational  amplifier  is 

necessary  to  provide  the  integrator  used  in  conjunction  with  the 

linear  accelerometer  (see  Fig.  1  for  the  block  diagram  of  the 

system).   Additionally,  an  accurate  tachometer  to  measure  the 

rotational  velocity  of  the  wheels  is  assumed  to  be  available. 

Thus,  all  elements  necessary  to  accurately  register  the  occurrence 

of  a  skid  are  assumed  to  be  available.   Our  attention  will  be 

directed  toward  the  actual  mechanism  by  which  the  triple  ramp 

control  is  generated,  and  by  which  the  pressure  applied  to  the 

brake  cylinder  is  modified. 

Figure  35  illustrates  the  means  by  which  it  is  proposed  to 
implement  the  triple  ramp  control.  An  explanation  of  Fig.  35 
follows . 

A  force  which  is  applied  to  the  brake  pedal  is  transmitted  to 
the  brake  drum  via  the  hydraulic  lines.   Until  a  skidding  situation 

occurs  V  is  held  equal  to  V    by  means  of  the  ramp  generator  and 

c  n        max 

limiter.  When  V  equals  V    ,  full  braking  action  is  applied  and 
c         max 

no  corrective  action  is  occurring. 

The  first  time  a  skid  occurs,  e  (error)  becomes  greater  than 
zero,  and  switches  (V)  and  (T)  throw  in  the  indicated  manner. 
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Fig.  35.  Proposed  implementation  of  tho  triple  ramp  control. 


Capacitor  C  discharges  rapidly  through  resistor  R~  (a  very  small 
resistance).  V  thus  drops  very  rapidly,  permitting  the  braking 
force  to  be  decreased  toward  zero;  this  step  is  an  exponential 
decay  with  a  very  small  time  constant,  but  it  approximates  the 
first  upward  ramp  of  F   (shown  with  a  slope  of  K„). 

When  the  error  is  reduced  to  zero  (skid  eliminated)  switches 
(l)  and  (z)  throw  once  more.   The  voltage  on  C   then  starts  to  climb 
with  a  slope  of  K, ,  thus  causing  the  braking  force  to  increase  once 
more.   During  this  first  part  of  the  cycle  a  portion  of  V   is  being 
monitored  by  C   (thru  the  tap  on  R,).  When  the  braking  force  has 
increased  sufficiently  to  cause  skidding  to  occur  once  again, 
switches  (V)    and  (2)  throw  once  more  and  V  drops  rapidly,  decreasing 
the  brake  force  as  before.  Additionally,  however,  switches  (\3)  and 

(2)  are  now  thrown,  and  capacitor  C  has  a  charge  proportional  to 
the  force  required  to  cause  skidding  (say  .9  x  V  to  cause  a  skid). 

When  V   falls  sufficiently  to  eliminate  the  skid,  e  =  0  again, 

and  the  ramp  generator  charges  C  upward  once  again.  When  V   is 

just  equal  to  the  charge  on  C   (.9  x  V   to  cause  a  skid,  for  example) 

V     becomes  greater  than  zero,  and  switch  (5)  throws  as  indicated, 
test         °  v— / 

This  puts  a  large  capacitance  (C.)  in  parallel  with  C„  and  causes 
the  ramp  slope  to  change  to  a  much  shallower  slope  (i.e.,  to  change 
from  slope  =  K,  to  slope  =  K  ).   V   is  thus  caused  to  approach  the 
value  that  creates  a  skid  much  more  slowly. 

When  (or  if)  a  skid  occurs  once  more,  the  entire  process 
described  above,  repeats  itself. 

A  sketch  of  the  force  acting  on  the  brake  cylinder  is  shown  in 
Fig.  36. 
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Fig.  36.   Force  acting  on  the  brake  cylinder  during  occurrence 
of  one  skid  pair,  if  using  Proposed  Implementation. 


Note  also  that  the  sensing  of  the  first,  third,  etc.,  skids 
or  second,  fourth,  etc.,  skids  can  be  accomplished  by  an  R  -  S 
flipflop  connected  to  the  output  of  a  comparator  which  would  be 
used  to  indicate  the  presence  or  absence  of  a  skid.   This  was 
done  on  the  analog  computer  simulation  with  very  satisfactory 
results . 

Figure  37  is  a  blow  up  of  the  pressure  controller  shown  in 
the  inset  of  Fig.  35.   This  is  a  standard  form  of  control  valve, 
and  is  not  original  with  the  Author. [l]  A  brief  explanation  of 
the  workings  of  a  standard  power  brake  system,  for  which  this 
particular  form  of  implementation  is  ideally  suited,  is  appropriate 
here. 

In  a  power  brake  unit,  when  the  brakes  are  applied,  a  valve 
to  the  atmosphere  is  opened  on  one  side  of  a  piston,  and  the  other 
side  is  connected  to  a  vacuum.   The  atmospheric  pressure  forces 
the  piston  closed,  and  causes  pressure  to  be  transmitted  to  the 
brake  cylinder  via  a  hydraulic  line. 
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Fig.  37.   Detail  of  pressure  controller  suitable  for  use 
with  a  power  brake  system. 


The  above  controller  will  cause  a  disruption  of  the  pressures 
on  the  power  cylinder  if  the  valve  is  rotated  from  the  position 
shown.   In  the  position  shown,  the  vacuum  line  is  blocked  off  thus 
permitting  full  atmospheric  pressure  to  the  piston;  any  other 
position  of  the  valve  would  lessen  the  pressure  on  the  piston, 
thus  reducing  the  braking  action. 
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9.   CONCLUSIONS. 

Of  the  three  major  types  of  controls  simulated,  the  triple 
ramp  control  proved  to  be  the  only  one  lacking  major  defects,  in 
that  it  reduced  the  degree  of  skid  to  a  very  low  level  while  still 
providing  a  very  short  stop  time.   Nevertheless,  the  double  ramp 
control  was  determined  to  provide  a  definite  improvement  over  the 
completely  uncontrolled  situation,  as  did  the  bang-bang  controls. 
However,  for  the  double  ramp  control  a  slightly  increased  stopping 
time  was  the  penalty  paid  for  obtaining  a  lower  degree  of  skid 
while  stopping  the  vehicle.   The  bang-bang  control  with  a  dead 
zone  proved  unsuitable  due  to  the  presence  of  a  slight  degree  of 
skid  during  the  entire  braking  path  of  the  vehicle.   Additionally, 
the  level  to  which  the  skid  could  be  reduced  would  be  severely 
limited  by  the  response  characteristics  of  the  brake  lines.   The 
bang-bang  control  with  a  time  delay,  while  still  superior  to  the 
uncorrected  situation,  appeared  to  be  considerably  inferior  to  any 
other  form  of  control  simulated. 

The  triple  ramp  control,  however,  appears  to  provide  very 
nearly  optimum  braking  conditions  in  that  the  degree  of  skid  and 
the  skid  time  are  both  minimized,  while  still  providing  a  considerably 
shorter  stop  time  than  would  be  obtained  from  any  other  form  of 
control.  A  trade  off  between  stop  time  and  adaptability  to  improving 
road  conditions  is  necessary  though,  and  the  proper  balance  between 
these  two  factors  would  be  determined  by  both  the  type  of  vehicle 
involved  and  the  anticipated  road  condition.   However,  greater 
emphasis  should  probably  be  placed  on  minimizing  skid  time  and  stop 
time  by  having  a  relatively  low  value  for  K  . 
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The  problem  of  possible  oscillation  introduced  by  the  triple 

ramp  control  should  never  arise,  since  it  is  essentially  a  negative 

[6] 
feedback  system  and  should  be  inherently  very  stable.   Extremely 

unusual  road  conditions  would  probably  be  necessary  for  any  oscil- 
lations to  occur;  however,  an  investigation  of  this  subject  has 
not  been  undertaken.   This  would  possibly  be  a  fruitful  area  for 
further  study. 

The  response  time  required  for  the  mechanical  components  of 
the  proposed  implementation  of  the  triple  ramp  control  (Fig.  35  and 
Fig.  37)  would  certainly  be  easily  attainable.   This  may  be  seen 
from  Fig.  25  (which  is  a  "typical"  response  curve),  where  examina- 
tion of  the  F  versus  time  curve  indicates  that  the  pressure  control- 
c  r 

ler  would  go  through  the  high  speed  cycle  at  the  rate  of  only  about 

4  c.p.s.   This  serves  to  give  an  order  of  magnitude  approximation 

to  the  response  characteristic  necessary. 

Throughout  this  entire  thesis,  F.   has  been  assumed  to  be  a 

in 

constant,  implying  a  panic  stop  situation  with  maximum  force  applied 

to  the  brake  pedal.   If,  during  the  course  of  slowing  the  vehicle, 

the  operator  decreased  the  pressure  applied  to  the  brake  pedal,  no 

effect  would  be  noticed  in  the  operation  of  the  automatic  control, 

as  lone  as  the  value  of  F.   were  still  sufficiently  large  to  cause 

in 

a  skid  (if  acting  by  itself).   If,  however,  F.   were  decreased  below 

the  value  needed  to  cause  a  skid  (if  acting  alone),  then  the  control 

would  permit  the  net  force  F  to  rise  to  the  new  value  of  F.   (either 

in 

at  the  rate  of  K,  or  K  ,  depending  upon  which  part  of  the  cycle  was 
acting  at  the  instant  of  decreasing  the  value  of  F.  )  and  this  value 
of  F  would  be  retained  until  a  skid  once  again  occurred.   If  the 
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applied  brake  force  F.   were  increased,  the  action  of  the  triple 

in  r 

ramp  control  would  be  unaffected,  as  the  control  would  still 
continue  to  seek  the  optimum  net  force  F. 

A  matter  of  possibly  more  concern  is  the  characteristic  of 
the  linear  accelerometer  and  of  the  integrator  used  to  indicate 
the  existence  of  a  skid.   The  proper  operation  of  this  skid  detect- 
ing scheme  is  basic  to  the  operation  of  any  skid  correcting  device 
proposed  in  this  thesis.   The  mounting  of  the  linear  accelerometer 
would  be  very  critical,  and  a  position  on  the  axle  housing,  in 
order  to  minimize  "nose  dive"  effects,  would  be  the  most  nearly 
ideal  location.   However,  no  serious  theoretical  limitations  on 
the  effectiveness  of  this  form  of  skid  detecting  scheme  present 
themselves,  and  this  would  appear  to  be  the  only  practical  means 
of  detecting  a  skid  prior  to  actual  wheel  lock  up. 

A  final  comment  concerning  the  analog  computer  on  which  the 
brake  system  and  all  controls  were  simulated  is  appropriate  at 
this  time.   The  COMCOR  Ci  5000  Analog  Computer  is  in  actuality  a 
form  of  hybrid  computer,  in  that  extensive  logic  circuitry  is  avail- 
able for  use  in  conjunction  with  standard  integrators  and  operational 
amplifiers.   Simulation  of  this  type  of  problem  would  be  difficult 
if  not  impossible  on  a  more  conventional  analog  computer  lacking 
the  logic  circuitry,  high  speed  switches  and  overall  versatility 
of  this  particular  computer.  A  simulation  of  the  type  described 
would  be  possible  on  a  digital  computer  utilizing  one  of  the 
languages  available  for  this  purpose,  such  as  DSL-360,  etc.   How- 
ever, in  order  to  obtain  the  extensive  data  needed  to  determine  the 
performance  characteristics  of  the  various  controls,  large  quantities 
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of  very  expensive  digital  computer  time  would  be  required.   On 
the  analog  computer,  however,  almost  two  hundred  different  runs 
were  made  consecutively  in  a  relatively  short  (and  inexpensive) 
time  period  by  merely  varying  potentiometer  settings  between  runs 
In  this  case,  the  accuracy  of  a  digital  computer  was  unnecessary, 
and  the  (hybrid)  analog  computer  proved  to  be  an  exceedingly 
valuable  design  aid.   This  thesis  thus  illustrates  graphically 
the  continuing  need  for  increasingly  more  sophisticated  analog 
computers,  such  as  the  COMCOR  Ci  5000. 
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APPENDIX  A 

Determination  of  the  maximum  possible  deceleration  of  a  vehicle. 

The  maximum  deceleration  possible  for  a  vehicle  is  determined 
by  the  coefficient  of  friction  between  the  tires  and  the  road 
surface.   Figure  Al  represents  an  idealized  vehicle  traveling 
in  the  plus  direction  with  a  velocity  X,  and  a  deceleration  X,  as 
indicated. 


Fig.  Al.   Idealized  vehicle. 

In  Figure  Al,  F  represents  the  force  of  friction  acting  on 
the  vehicle,  and  is  the  only  external  force  causing  the  vehicle 
to  decelerate.   The  weight,  mg,  is  the  normal  force  acting  on 
the  ground,  so  the  value  of  Ff  is  given  by 

Ff  =  Mmg  }  (Al) 

where  p  is  the  coefficient  of  friction  between  the  wheel  and  the 
road  surface.   The  coefficient  of  friction  is  a  variable,  which 
has  its  greatest  value  (p  )  just  before  a  skid  occurs,  and  then 
drops  abruptly  to  a  lower  value  |i  when  a  skid  begins. 
Thus  the  maximum  value  of  F   is 

Ff  (max)  "  ^  .  <A2> 
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Note  that  |Jb  corresponds  to  the  coefficient  of  kinetic  friction, 
whereas  \i     corresponds  to  static  friction. 

Summing  the  forces  in  the  X  direction,  we  obtain 

mX  =  -  Ff  .  (A3) 

Substituting  Eq.  (A2)  in  Eq .  (A3)  and  solving  for  X.    . ,  we 

(max) 

obtain  for  the  maximum  possible  deceleration 


"  X(max)  "  V  •  <A4> 
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APPENDIX  B 


Explanation  of  analoc  computer  ggbgla  used. 

The  analog  computer  elements  used  in  the  simulations  described 
in  this  thesis  are  indicated  below. 


SYMBOL 


out 


X*- 


X, 


Xoot 


X—  x 


out 


<& 


Xout 


oXTe5T>  K 


yr-~* Xout 


4 o 


Xt«t>K 


~Y- 


■nut 


c 
c 


D£3CRIPnOM 
Sign  changing  summer  amplifier. 

Xout^    ~(xi  "HX£   +  I0X3) 

Sign  changing  integrator. 


r  t« 


Xout 


f  (x,-HOX£)  c/t    +X< 


Claas  HI  Multiplier. 

Xout  =  X  i  X  g 

Potentiometer . 

Xout*=  K,x,  ,  K,<  I 

Double  pole,  double  throw  relay. 

Xou-t  =  X  |  ;  "for   XTeST  <C  K 

Xout  =  X£  3  for  Xtest  >  K 
High  speed  digital/ analog  switch. 

Xout  =  O    ;-for     XTetT<K 
Xout=    X,   j  -For     XTetT>  K 

Comparator. 

C  *  logical  1  and  C  «  logical  0, 

when  %i  +  x2>0. 

C  ■  logical  0  and  C  *  logical  1, 

when  x^  ♦  X2<0. 

Comparators  are  used  to  throw  both 

the  DPDT  relays  and  the  D/A  switches. 
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-Y- 


i        Point 
1 Xouf 


*^H  - 


TF5T(>K,  TFST£>Ke 


Xo«t 


Tio  point. 

Used  in  conjunction  with  two  or  «oro 

D/A  switches  which  are  thrown  so  that 

only  on*  switch  is  ON  at  a  tiae. 

Xou-t  =  X)  or  x^ 

Switching  accomplished  by  HAND  logic, 

utilising  only  one  D/A  switch. 

Xout  =  X,  j     when 

CTEST,>K,)   anol     (TEST£>ka). 

A  Iso  ,  Xou-t  =  Oj  when 

(TEST,  <K,)   or     (TEST2<K2). 
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